function [D,h_D] = get_circle2(contours,i,M,dr,kh)
% difference from get_circle is that dr is percent of minimal radius to add
C        = contours{i};

xc       = mean(C(1,:));
yc       = mean(C(4,:));
R        = max(sqrt((C(1,:) - xc).^2 + (C(4,:) - yc).^2 ));

theta    = linspace(0,2*pi*(1-1/M),M);
D        = zeros(2,M);
RR       = R*(1 + dr);
D(1,:)   = xc + RR*cos(theta);
D(2,:)   = yc + RR*sin(theta);
h_D      = 2*pi*RR / M;

% fprintf(1,'Outer contour D1 \n')
% fprintf(1,'  center : (%f ,%f)\n',xc,yc)
% fprintf(1,'  radius : %f\n',RR)

return